COMPUTERS 


7524 S.W. MACADAM AVE. + PORTLAND 19, OREGON ‘MARCH 


REPLACES AUGUST 1959 


AN ILLUSTRATIVE PROBLEM FOR THE ESIAC COMPUTER 


DESIGN OF TACHOMETER STABILIZATION 


In a positional control system for a heavy inertia load, the load is driven by a 
dc gear motor which receives its power from an Amplidyne power supply. A 
synchro transmitter is connected to the load and a synchro control trans- 
former to the input control shaft. The output of the synchro control trans- 
former is synchronously rectified in a phase sensitive demodulator consisting 
of a push-pull pair of 6L6 tubes using the carrier frequency for the plate 
supply. The output of the demodulator supplies the Amplidyne field. 


The performance of the system requires that the open loop gain coefficient be 
adjusted to a value of 9000. The ESIAC is used to show that at this gain set- 
ting, the basic system as described above is unstable. The system will there- 
fore be modified by the addition of tachometer feedback. The block diagram 
and equations for this system are shown in Figure l. 
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The purpose of this analysis is to determine the optimum tachometer 
coefficient. The method used will be to determine first the locations of the 
closed loop poles of the system without tachometer feedback (when K = 0), 
then to determine the change in the locations of the poles as a function of K_, 
selecting the value of K, which will yield the maximum damping of the systém. 
The locations of the polés will be determined by root locus techniques as 
plotted by the ESIAC. 


The overall transfer function of the system is 


6 
Se : 
eo. s(s + 6) (s + 12) 
1 Kis + 1 +9000 
This equation can also be written 
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in which s., s, and s, are the roots of the equation 


3° "4 5 
s(s + 6) (s +12) _ 
1 + "9000 a 


or s(s + 6) (s + 12) = - 9000 


This equation corresponds to the ESIAC equation in the form 


Ny i a 
F = Ks Tf (s - s,) 
in which F = 9000 2180 
K = +1 
no = + 1 
8) = 6£180°,n, = +1 
s> =12 2180 »n,=ti 


The steps in solving this equation on the ESIAC are as follows: 


1. Set the controls on the right hand panel so that the equation reads 


2. Set IFI = 9000 
3. Set ZF = 180° 


4. SetK=1 
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5. Plot points s = 6Z180° =s, ands = 122180° =s. ona sheet of 
ESIAC graph paper. (This graph paper has a logarithmic |s! scale 
and a linear Zs scale corresponding to the ESIAC scales.) Place this 
graph sheet on the plotting board of the ESIAC. 


6. Move the plotting stylus to point s, on the graph paper. Place one 
probe at the index point at the center of the phase sheet and another 
probe at the index point at the center of the magnitude sheet, using 
probe clamps of value 1 screwed into holes marked +. 


7. Repeat step 6 for point S5. 
8. Move the plotting stylus in a small circle around point s, to finda 
point where the ZF meter indicates a null, then move the stylus 

away from the point along a path where the ZF null is maintained, 


marking any point where the IF| meter also indicates a null. Each 
such point is a root of the equation. 


9. Repeat step 9 starting from the point S5. 


10. Repeat step 9 starting from each point where an ZF null can be found 
along the left and right borders of the graph paper. 


The plot resulting from the above procedure, with the curves along which 
4F is at a null drawn in as solid lines, is shown in Figure 2. 
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Figure 2 
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There are three roots to the equation, located at 


s, = 18276" 
8, = 18 2-76" 
sp = 27 Z180° 


These roots are the poles of the system without tachometer feedback, i.e., 
when K,=0. It is. seen that roots s, ands, represent an unstable response 
: no iv 
(an increasing oscillation), since their angies are less than 90°. Modifi- 

cation of the system is necessary in order to make it operable. 


Using the results of the above calculation, the transfer function of the 
system with tachometer feedback can be written as 


a5 _ l 
0; Ks + (1 - ————) (1 - ——~—_5 (1 - —+—__ 
18276° 18Z - 76° 272180° 
This equation can also be written 
So. 
6. S Ss S 
i (Loe) thes) (hee 
6 7 8 


in which Sp: S. and S. are the roots of the equation 
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This equation corresponds to the ESIAC equation in the form 


el ra 
F = Ks II (1 = Ps ) 
i 
in which 
F = 121809 
K = Parameter to be determined = KY 
n =+#+1 
fe) 
Oo. _— 
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The parameter K, affects only IFIl, not ZF. The method of solution will there- 
fore be to plot the loci of values of s which yield the specified value of ZF, 
then to select the points on these curves where the complex pair of roots 
correspond to a response with maximum damping (the points where Zs is 
nearest 180°) and at these points to determine the value of K, yielding the 
specified value of IFI. 


The steps in solving this part of the problem on the ESIAC are as follows: 


1. Set the controls on the right hand panel so that the equation reads 


l 


F=a¢ Ks" m(1-2) 
: 


2. Set IFIl=1 
3. Set ZF = 180° 


4, Lay anew sheet of ESIAC graph paper over the sheet previously 
plotted and mark the locations of the roots of the previous equation 
$3254 and Spe Place this graph sheet on the plotting board. 

5. Move the plotting stylus to point s,. Place one probe at the index 
point at the center of the phase shéet and another probe at the index 
point at the center of the magnitude sheet, using probe clamps of 
value 1 screwed into holes marked -. 


6. Repeat step 5 for point Sys 

7. Repeat step 5 for point Spe 

8. Move the plotting stylus in a small circle around point s, to find a 
point where the ZF meter indicates a null, then move thé stylus 
away from the point along a path where the ZF null is maintained, 
marking a sufficient number of points along this path to establish the 
curve. 


9. Repeat step 9 starting from point S 4° 

10. Repeat step 9 starting from point Spe 

11. Repeat step 9 starting from each point where an ZF null can be found 
along the left and right border of the graph paper. 


12, Locate and mark the points where the curves starting from s, and s 
approach nearest the line along which Zs = 180°; this position of the 
roots yields maximum damping. These roots are located at 


4 


s 


41 299° -6 + j40 
41Z-99° = -6 - j40 


ay 
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13. Adjust the K control so that a IFl null is indicated. Read the value of 
K,(0.18). This is the value of the tachometer feedback coefficient 
which is the solution to the problem. 


14. To locate the other root of the equation for use in analysis of the 
system performance, move the stylus along the path starting from 


point s, until the point is found where a IF! nullis indicated. This 
root will be found at 
Sg = 5.42180 = -5.4 


The plot resulting from the above procedure is shown in Figure 3. 
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Figure 3 


The overall transfer function of the system is 


0. Ks+1+S8tonet (i < Ss - (1 - S (1 - Ss : 
41 299 41 2-99 5.42180 


As an indication of the performance with this value of K,, the frequency 
response curves (gain and phase shift) are easily plotted on the ESIAC. The 
overall transfer function of the system when s = jw is 


9 
OL I 
0. jw _ jw _ __ jw 
: (1 - g357990) (1 - GxZ 990) (1 - EGTA) 
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This equation corresponds to the ESIAC equation in the form 


8 
in which IFl = |g-|= gain 
i ; 
5 
ZF = £—2= phase shift 
Oi 
= +1 
n = QO 
1°) 
s, = 41499", n,=- 1 
6 ae ae 
s7_ = 412-99", n =- 1 
Oo 
Sg = 5.424180 , Ny =- l 


The steps in plotting the gain and phase shift of the system on the ESIAC are 


s =jw =0Z90° 


as follows: 


1, 


3. 


Set the controls on the right hand panel so that the equation reads 
: \ 


F= +Ks° I (l- =) 
i 


SetK=1. 


Place on the plotting board the graph sheet on which s,, s_ ands 
> 6° 7 8 
have been plotted 


Move the plotting stylus to point s,. Place one probe at the index 
point at the center of the phase sheet and another probe at the index 


point at the center of the magnitude sheet, using probe clamps of 


value 1 screwed into holes marked -. 


| Repeat step 4 for point Sa. 


Repeat step 4 for point Sg- 


Remove the graph sheet locating points Sy, $7 and Sq. 


Set Zs = 90°, 
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9. To plot gain, couple the plotting mechanism for IFI versus Isl . 
Place on the plotting board a sheet of ESIAC graph paper for gain 
plots. Move the plotting stylus in the Isl direction to any selected 
value of Isl = and then in the IF! direction until the [Fl meter 
indicates a null. Mark the point. Repeat at as many values of w 
as it is desired to plot for the curve. The resulting graph is shown 
in Figure 4. 
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10. To plot phase shift, couple the plotting mechanism for ZF versus Isl . 
Place on the plotting board a sheet of ESIAC graph paper for phase 
plots. Move the plotting stylus in the Isl direction to any selected 
value of Isl] =a and then in the ZF direction until the ZF meter indi- 
cates a null. Mark the point. Repeat at as many values of w as it is 
desired to plot for the curve. The resulting graph is shown in Figure 5. 
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The time domain response for any specified input, 0, may be found by 
taking the inverse Laplace transform of 0.. The ESIAC can be used to find 
the residues at the poles of 00 and then the time domain response may be 


written by inspection. | 


For example, assuming that the input is a unit step 


rewrite the transfer function 


wa 7 9000 
oO = te ea ae oe te es) 
0. (s S¢ S S27 S 5 9) 
then 
9 = 9000: 
re) ‘gua: S- 8. 5 - 59) § 
or 
A A A A 
8, = soe teeta te 
6 i a7 Pg 
in which Ag: Av, Ag and Ay are the residues at Ses Sa: 89 and 0 


These are found on the ESIAC to be 


A, = 0.07280" at s, = - 6+ j40 
A, = 0.072-80° at 8, = + 6- j40 
o = on 
Ag = 12180 at s, = - 5.4 
and by inspection 
AS = | at s = Q 


Taking the inverse Laplace transform term by term the output as a function 


of time is 


bt * 360° -5. 4t 


= O 2 
8 =f 6, = 0.14 e cos (80 fe AOr ae) ace + 1 
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This problem has illustrated,several applications of the ESIAC which may be 
summarized as follows: 


1, 


FACTORING: The fixed portion of the denominator (the part excluding 
the K, term) was factored by means of the root locus technique on the 


ESIAC. 


ROOT LOCUS PLOTTING: The locus of the closed loop poles of the 
system as a function of K, was made on the ESIAC, and the optimum 
pole locations and corresponding value of K, determined. 


FREQUENCY RESPONSE PLOTTING: Plots of gain and phase as a 
function of frequency (Bode plots) for the system were made on the 
ESIAC. 


TRANSIENT RESPONSE CALCULATION: The Laplace transform of 
the system output for a specified input was expanded as the sum of 
elementary terms and the coefficient of each term (the residue) was 
determined on the ESIAC. With this information the output as a 
function of time was written by inspection. 
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